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ABSTRACT 

A simple, semi- analytic representation is developed for nuclear burning in Type la supernovae in the 
special case where turbulent eddies completely disrupt the flame. The speed and width of the "dis- 
tributed" flame front are derived. For the conditions considered, the burning front can be considered 
as a turbulent flame brush composed of corrugated sheets of well-mixed flames. These flames are 
assumed to have a quasi-steady-state structure similar to the laminar flame structure, but controlled 
by turbulent diffusion. Detonations cannot appear in the system as long as distributed flames are 
still quasi-steady-state, but this condition is violated when the distributed flame width becomes com- 
parable to the size of largest turbulent eddies. When this happens, a transition to detonation may 
occur. For current best estimates of the turbulent energy, the most likely density for the transition to 
detonation is in the range 0.5 — 1.5 x 10 7 g cm~ 3 . 

Subject headings: supernovae: general; hydrodynamics, shock waves, turbulence 



1. INTRODUCTION 

One of the greatest uncertainties in how a Chan- 
drasekhar mass white dwarf explodes as a Type la su- 
pernova is whether and how an initially subsonic burn- 
ing front, a deflagration, makes a transition to a super- 
sonic detonation. A related question is the characteris- 
tics of nuclear burning in a medium where turbulence 
has become so strong that hot ash and cold fuel can 
be co-mingled before burning. The conditions for the 
latter, known as "burning in the distributed regime", 
have long been known to both the combustion and 
astronomical communities O Khokhlov. Or an. fc Wheeler 
199ft INiemever fc Kerstein |199ft iNiemever fc Woosley 
199ft lPeterdl2000l iPope Ill987f ) . When the laminar flame 



grows thick enough and the turbulent intensity great 
enough, the "Gibson" length, that length scale which 
turns over due to turbulent eddies as fast as a laminar 
flame crosses it, becomes smaller than the flame thick- 
ness itself. That is, 
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Approximate conditions for entering the distributed 
regime ha ve been given for an explod ing white dwarf in 
Fig. 2 of INiemever fc Wooslevl (|1997| ) and those condi- 
tions define the applicability of this paper. For those con- 
ditions, icib is also much greater than the Kolmogorov 
length, L Kol = LRe~ 3 ^ ~ 10~ 4 - 1CT 5 cm, where the 
turbulence is dissipated. Here L is the integral length 
scale of turbulence (~ 10 6 cm) and Re is the Reynolds 
number (~ 10 ). The inequality LkoI -C ^Gib <C <5iam is 
thus satisfied, where <5i am is the thickness of the laminar 
flame (typically ~cm). 

Deep in the distributed regime, turbulence is more ef- 
fective at transporting both heat and composition, even 
on scales as small as the laminar flame width, than con- 
duction and diffusion (the c onductive and r adiative opac- 
ities here are comparable; iTimm es 2000). As a con- 
sequence, the concept of a laminar flame, one whose 
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width and speed are determined by the near equality 
of bu r ning and diffusion time sc ales ([Landau fc Lifshita 
119591 : ITimmes fc Wooslevl Il992h . breaks down. Fuel 
and hot ash are co-mingled and the definition of th e 
width of the must be modified (jLisewski et al J 120001) . 
If that width grows large enough and the burning 
rapid enough, a transition to detonation can occur 
dKhokhlov. Oran. fc Whee^[l997l: iLisewski et alJl200Ct 
INiemever fc Wooslevlll997| ). Such a transition is impos- 
sible within the extent of a laminar flame; it can only 
occur in a turbulently stirred one. 

Here the physics of that transition is explored. Kol- 
mogorov scaling is assumed throughout. Expressions 
are derived for the turbulent flame speed and its width 
(§[2). A hypothetical "steady state" is posited, in which 
fuel burns at a rate balancing turbulent mixing. The 
width of the flame grows as the density declines be- 
cause the temperature of the ash, to which the burn- 
ing rate in the mixture is very sensitive, is lower there. 
The turbulent flame moves at a rate given by the length 
scale at which turbulent diffusion matches burning and, 
as the width of the burning mixture becomes greater, 
so too does its speed. Thus, in supernovae, the dis- 
tributed burning flame moves faster as the density de- 
creases - the opposite of what happens for a laminar 
flame. On larger scales, turbulence also folds these (tur- 
bulently broadened) flames and, just as in the laminar 
case, there is an overall "flame brush" (jDa mkohler l l 19401 : 
iHillebrandt fc Niemeverlll999l ) whose motion determines 
the total rate of burning. 

As the density continues to decline, however, the tur- 
bulently mixed flame grows ever broader, eventually ap- 
proaching the size of the integral length scale of the 
turbulence. Qualitatively, this is the largest scale at 
which the anisotropic shear and instabilities introduced 
by floatation produce an nearly isotropic cascade of tur- 
bulence with constant energy density. Technically, it is 
the distance scale beyond which the self correlation of 
the velocity components vanishes. For a typical Type 
la explosion, plumes of size ^100 km float at speeds 
in excess of 1000 km s _1 . Empirically, from numerical 
simulations, the size of the isotropically stirred region is 
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~10 km, though certainly variations of a factor of several 
around this are allowed. Also important is the velocity 
at the integral length scale, ul, which specifies the en- 
ergy density in the turbulence. This is typically some 
fraction, ~10%, of the floatation speed, but again, large 
variations are expected, up to the floatation speed itself. 

When the width of the turbulently mixed flame be- 
comes comparable to the integral length scale of the tur- 
bulence itself, the largest turbulent speeds have thus be- 
come an appre ciable f ractio n of the sound speed. Re- 
cent studies bv lRbpke] (|2007l ) show the highest turbulent 
speeds can approach 1000 km s _1 at densities ~10 7 g 
" 3 ; the sound speed there is about 4000 - 5000 km 
So long as the assumed steady state solution per- 
sists, supersonic burning remains, by definition, impossi- 
ble for subsonic turbulence. However, once there is only 
a single flame or two in the integral length scale, the 
steady state assumption certainly breaks down. A large 
eddy and its accompanying cascade stirs the mixture, 
but then burning goes on at a non-steady rate before 
another eddy happens in the same region. In § 14.21 and 
§ l4.3i it is shown that this sets the stage for a detonation. 
The conditions are restrictive and require turbulent en- 
ergies corresponding to close to 1000 km s _1 on a length 
scale of 10 km and a density between 0.5 and 1.5 x 10 7 
g cm -3 . This is a much more restrictive condition than 
just "entering the distributed burning regime" . 

This work bears some simil arity to 
those of |Khokhlov. Oran. fc Wheelerl 
iNiemever fc Wooslevl (|1997t ). and iLisewski et ail 
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but considers more carefully both the conditions in the 
mixed region and the need for a turbulent flame that is 
already moving at a fraction of the sound speed before 
a transition to detonation can occur. The conditions 
required are thus more precisely determined and much 
more constrained. The necessity of a separation between 
the carbon burning flame and the oxygen burning flame 
and the breaking of the steady state assumption at 
the integral length scale are emph asized. S o me of the 
conclusions are similar to those of iKers tcin (200l]) for 
low Prandtl number flames, but are more extensively 
discussed in the astrophysical context. 

2. A SIMPLE MODEL FOR BURNING IN THE 
DISTRIBUTED REGIME 

The basic parameters of any flame, its width and speed, 
can be estimated from the rates for fuel consumption and 
fuel advection into the burning region. Although turbu- 
lence is stochastic and the distribution of heat in a stirred 
region is not nearly so smooth as when huge numbers of 
electrons and photons participate in conduction and dif- 
fusion, there is still, on the average, something like a 
steady state. In that steady state, the rate at which fuel 
(carbon) is brought into the burning region by turbu- 
lent eddies balances the rate of consumption within that 
region. The whole mixed-up, burning ensemble moves 
through the fuel with a typical speed which is essentially 
the size of the region divided by its turbulent turnover 
time. That is, 
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where nyi = pN aY\2 is the number density of carbon nu- 
clei as a function of location, p is the density, Na, Avo- 
gadro's number, Y±2, the mass fraction of carbon, X12, 



divided by 12, Vt, the velocity of the burning turbulent 
front normal to the area, A, and V , the volume bounded 
by that area. Here we consider a one-dimensional flame 
in plane geometry. For the low densities of interest, reac- 
tions beyond carbon burning, i.e., oxygen burning, occur 
so far behind the carbon burning flame as to be negligi- 
ble on the scale of the problem. Let the thickness of the 
mixed region be A, then 
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The left hand side gives the rate of carbon destruction 
(in gm cm -2 s _1 ) by nuclear reactions, while the right 
hand side is the rate of advection into the burning region 
by turbulent eddies with characteristic length scale, A, 
and speed Vt- Here, pf uo i is the density in the unmixed 
fuel and, Y° 2 , the carbon abundance there. The burning 
rate is 

dY 

= -2 P y 2 i? 12 , 12 (p,T), (4) 

where i?i2,i2 is the rate factor for the carbon fusion reac- 
tion. Because of the temperature sensitivity of this rate, 
about T 20 , most of the carbon consumption goes on in 
a narrow region w here the mass fr action is low and the 
temperature high (|Bell et al.ll2004h . closer to the hot ash 
than to the cold fuel. The rapid rate of burning there 
has to balance, on the average, what is advected into the 
larger region. In order to obtain vt, it is thus necessary 
to specify p(l'), T(l'), Y 12 (l'), and v T - 

For vt, it is appropriate to take the velocity of the tur- 
bulent eddy with length scale A. Assuming Kolmogorov 
scaling and a turbulent energy input on the integral 
length scale, L, corresponding to velocity ul, 
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where, from typical numerical simulati ons, Mr, is in the 
range 10 7 - 10 8 cm s" 1 for L = 10 6 cm (|R5pkdl2007l) . 

All speeds are very subsonic, so to good approximation, 
the pressure in the fuel, ash, and the mixture is constant. 
Provided that mixing is faster than burning and conduc- 
tion, the temperature and density in the mixture are thus 
uniquely defined by the fractions of ash and fuel that 
are mixed and the compositions of each. In fact, some 
burning does occur during the mixing and this affects 
the temperature distribution, but to first order, the tem- 
perature obtained by burning to a certain carbon mass 
fraction is the same as that obtained by mixing cold fuel 
with ash of higher energy and lower carbon abundance to 
obtain that same final mass fraction. This is not precise 
because carbon burns to different products at different 
temperatures and so the energy released is not a linear 
function of carbon consumed, but the difference is not 
large. 

Next we seek a description of how temperature and 
density vary in the mixed region. This requires an 
ansatz for how the carbon mass fraction varies. To il- 
lustrate the procedure, assume an initial composition of 
50% C and 50% O at a density of 1.0 x 10 7 g cm" 3 
and temperature 6.0 x 10 s K. This is a typical tem- 
perature in the outer parts of the white dwarf when it 
runs away, but the answer will not depend on the ex- 
act value because the pressure is not very sensitive to 
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the temperature and the internal energy of the ash is 
much higher than that of the fuel. The pressure in this 
fuel is 9.046 x 10 23 dyne cm~ 2 and its internal energy is 
e = 1.706 x 10 17 erg g _1 (here and throughout the pa- 
per we employ the Helm holtz equation-of-state routine of 
iTimmes fc Swestvll2000( ) . Now mix this fuel with a small 
amount of ash so that the temperature rises a small in- 
crement, ST, and the carbon fraction goes down. For 
this new state, (T,P), iterate on the density and internal 
energy until a solution is found with the same pressure 
as before. The new density is lower and its formation re- 
quired expansion. Its new energy is the heat brought in 
by the mixing minus the energy lost to PdV work. That 
is 
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e(T a + ST) - e(T ) - P(H). 



Here dp is inherently negative so the pressure term is 
positive. The composition of the mixture will also have 
changed by an amount that depends upon the composi- 
tion of the ash. Here we adopt the ash composition given 
by following isobaric burning to completion off-line us- 
ing a small 7 isotope network. For the conditions given 
above, a typical ash composition will be 57% O, 16% Mg, 
26% Si and 1% S. At a density of 3 x 10 7 g cm" 3 , the 
composition would have been slightly different: 57% O, 
8% Mg, 33% Si, and 2% S, but the Q-value is not very 
sensitive to the difference. The change in nuclear binding 
energy between the ash and fuel (50% each 12 C and 16 0) 
is thus Q = 3.10 x 10 17 erg g _1 . Substantially different 
numbers characterize an initially carbon-rich composi- 
tion. For a fuel that is 75% carbon and 25% oxygen, the 
ash composition at 10 7 g cm' 3 is 38% O, 15% Mg, 45% 
Si, and 2% S, implying Q = 4.70 x 10 17 erg g" 1 . 
The fraction of carbon in the mixture is 



Y 12 (T + 5T)=Y 12 (fuel)(l - % 
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since, by definition, Yi2(ash) — 0. The other composition 
variables are similarly interpolated between their initial 
(fuel) and final (ash) values based upon the change in 
energy. Using this new composition, the density is again 
iterated to find the isobaric state appropriate to the new 
temperature and self-consistent mixed composition. The 
process is continued for about 1000 steps until the carbon 
abundance is zero. The outcome of this calculation is a 
set of temperatures, densities, and compositions for the 
mixture consistent with the pressure in the fuel. 

To reach closure, it remains to specify how the carbon 
abundance varies within the mixed region. In reality, 
the carbon mass fraction will be heterogeneous, reflect- 
ing the operation of numerous eddies on all scales and 
the large Lewis number. Stirring is more effective on 
small scales so the most natural distribution would be a 
"nois y" staircase fu nction, or even a homogeneous mix- 
ture (|Kersteinll2001[ ). We return to this picture in § 14.21 
For now, however, a simple approximation is made that 
the carbon abundance is distributed linearly within the 
stirred flame. That is, for < I' /A < 1, 



(8) 



Such a linear approximation is consistent with multi- 
dimensional simulations so far (Fig. 31 of Bel l et all 



l2004h and gives equations that are easy to manipulate 
and understand. 

Given Yi 2 (l'), p(Yi 2 ) and T(Yi 2 ), L, and ul, one is now 
equipped to solve eq. ([3]) for a unique value of A. Some 
results are given in Fig.[T]and Table 1, the latter showing 
a dramatic dependence of the burning front width and 
speed on the density and turbulent energy. For the lowest 
densities and highest turbulent energies considered, the 
mixed region becomes comparable to the integral length 
scale, 10 km, and the burning can approach a fraction of 
the sound speed. The numbers which give As 10 km in 
Table 1 are not physical unless eq. (5) and the assump- 
tion of homogeneous, isotropic turbulence can be extrap- 
olated to these larger scales, and the velocity continues to 
increase above u^. In those cases where ul is already 10 s 
cm s -1 , that is doubtful. The flame widths and speeds 
are also smaller for carbon-rich mixtures. Burning more 
carbon raises the temperature of the ash and mixture 
and makes the burning region smaller. 

It is important that the flame has separated into two 
components. If one added the energy generation from 
oxygen burning, the temperatures would be higher and 
the width of the burning region much smaller. Since 
as we shall see, only the largest values for A in Table 
1 imply a possible transition to detonation, a necessary 
condition for a delayed detonation transition is that the 
oxygen and carbon burning flames have split and are 
widely separated. 

3. APPROXIMATIONS TO THE TURBULENT FLAME 
SPEED 

The speed of a laminar flame is proportional to the 
square root of the heat diffusion coefficient. In the dis- 
tributed regime one expects a similar relation with the 
turbulent diffusion coefficie nt, Dtmh ~ v\\ substitutin g 
for the radiative one (e.g., iRopke fc Hillebrandtl [20051. 
Here v\ is the turbulent eddy speed on the scale of the 
flame width, A, and hence D tU rh ~ A^u^L -1 / 3 . The 
width of the flame is given by equating the nuclear and 
diffusion times 

A 2 LV3 A 2/3 



and hence 
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This relatio n is well known in the chemical combustion 
community (|Kersteinll200H lPeterslfl999l . |2000[) . Here r nuc 
is the average nuclear burning time in the region defined 
by eq. (2). For a given density and turbulent energy, 

this suggests a scaling A oc u 3 ^ 2 which agrees with the 
values in Table 1. Such a scaling is also expected because 
the left hand side of eq. (|3|) is proportional to A while 
the right hand side depends on A 1 / 3 ul. Equation pop 
also states that the turbulent flame speed, A/r nuc , is the 
square root of the energy dissipated by the turbulent cas- 
cade, u\/ ' L in a nuclear time scale. For a given turbulent 
energy, to get the flame to move faster one must increase 
the nuclear time scale, i.e., slow the burning. 

The nuclear time scale is very sensitive to the temper- 
ature in the mixed region which is highly variable, but 
near 3 x 10 9 K, it is approximately ([Woosley et al.ll2004h 

~ CpT nu 

T„uc » , (11) 
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where 5 nU c °c pX\ 2 T n . The heat capacity, due to a com- 
bination of semi-degenerate electrons and radiation, in- 
creases very roughly as T 2 while, for barrier penetration, 
n in the temperature range near 2 to 3 billion K is 

n = 28.05 T~ 1/3 - - « 20. (12) 

o 

Thus r nuc is roughly proportional to p~ x A^ 2 T~ 17 and, 
since the temperature in the burning region is propor- 
tional to T as h, the flame width, 

A oc uTL-^p^T^X^ (13) 

~ 40 ( Zkh \-25.Bf u L \3/2 

V2.79xlO u V 10* cm s- 1 / 

(10kmU/2 ( Pfud )-3/2 ( 05_)3 k 
V L / V 10' g cm J / V A12 ' ' 

where the proportionality has been normalized to the 
numerical results in Table 1 for the fiducial values. This 
is an overall good fit for other densities and compositions 
in the range of interest. 

The turbulent speed, vt, is given either by eq. $5$ eval- 
uated for the flame width derived above, or by 

~ L 1 / 2 

~ 7Qn ( T aah \-8.5f ul ^3/2 

~ ,inJ V2.79xl0« K / V 10 s cm s- 1 > 

In addition, vt should not be greater than u^. Note the 
strong reciprocal dependence on the temperature of the 
ash and hence on the density and energy yield of the 
burning. It is this dependence that is chiefly responsi- 
ble for the spreading and acceleration of the turbulent 
flame at low density. If oxygen burned as well as carbon, 
T as h would be much greater and the flame speeds and 
widths would be drastically reduced. The fast speeds and 
broad widths derived here rely on the oxygen flame lag- 
ging far behind the carbon flame, i.e., outside the mixed 
region. For the low densities we consider, this is the case. 
Note also that vt is enormously greater than the laminar 
flame speed, wh ich at 10 7 g cm -3 is only 3000 cm s _1 
(|Bell et al.H200l . 

The quantity T as h can be computed off line 
as a function of initial density and composition 
(T immes fc Wooslcvi I1992D . Approximate values ob- 
tained using a small 7 isotope network are given in Table 
2. The values in the table can be approximately fit by 
an expression of the form 

/ v \ °- 25 

T ash « 2.79 x 10 9 - ^ 3 ) . (15) 

\5 x 10 b gem 7 

Substituting this in eq. (fT4|) and eq. (fT5|) . one has 

A « 4.9 ( 1Q8 Z ^0 3/2(10^)1/2 (16) 

( ln 7 Pfucl -, r 7 - 9 (^) 9 - 4 km, 

V 10' g cm d / V A 12 ' ' 

and 

790 (Tf^^r) 3 / 2 (17) 
( _ 7 ^ ) -2.6 ( 0^ ) 3.1 ( 10^ m)1 /2 kmg -l. 

For a given turbulent energy and carbon fraction, this 
implies a flame speed that scales roughly as p~ 26 



4. THE TRANSITION FROM A DEFLAGRATION TO A 
DETONATION 

4.1. Spontaneous Detonation 

One of the simplest ways a transition to detonation 
could happen in an exploding white dwarf, which is in- 
cluded here only because it seems to have been over- 
looked, is if frictional heating - the dissipation of the 
turbulent energy on the Kolmogorov scale - heats a re- 
gion of fuel to the flash point. 

Consider a region of fuel close to a large rising ele- 
ment of ash. The rise of a burning plume injects tur- 
bulent energy at some characteristic length scale. The 
velocities are initially anisotropic, but after the energy 
cascades down approximately one decade in length scale, 
that energy resides in isotropic Kolmogorov turbulence 
(jZingale et al.l l2005). Typical turbulent speeds at length 
scales of ~10 km where the Kolmogorov cascade might 
begin are, at reasonably late times in the explosion, in 
the ra nge 1 - 10 xlO 7 cm s" 1 (|R5pkdl2007tlSchmidt et all 
l2006lh 

This energy cascades downwards to the Kolmogorov 
length, 10~ 4 - 10 -5 cm, where it dissipates as heat. The 
amount of heat dissipated is u\jL erg g s , i.e., the con- 
served quantity in Kolmogorov turbulence. For speeds 
10 7 - 10 s cm s _1 on a length scale of 10 km, that corre- 
sponds to 10 15 - 10 18 erg g s -1 . Most of th is dissipation 
occurs inside the ash (|Schmidt et al.|[2006l ). in part be- 
cause the flame spreads at a speed comparable to that 
of largest eddies. However, there may be small regions 
near the floating ash, unresolved in current studies, per- 
haps within the Kelvin-Hclmholtz rolls that bound the 
rising plumes or in the wakes of detached bubbles, where 
a locally large concentration of turbulent energy is dissi- 
pated in the fuel. The dissipation might be particularly 
large in vortex tubes shed by the rising plumes. 

Because of its low heat capacity, s2 x 10 17 erg g _1 
is necessary to raise the fuel temperature to the point 
where it will burn supersonically. If the explosive burn- 
ing region is larger than a critical mass, a detonation 
will o ccur (jDursi fc Timmed 120061 : iNiemever fc Wooslevi 
Il997f h For the highest turbulent energies considered, this 
would only take about 0.2 s of uninterrupted dissipation, 
significantly less than the expansion time of the star. For 
regions in close proximity to an active flame, the friction- 
ally heated fuel will probably be burned before it can run 
away. However, in the trailing wake of rising bubbles, 
there might be time for viscous dissipation in unburned 
fuel to ignite new burning. If the temperature gradient 
in this new region is sufficiently small, the ignition could 
have a supersonic phase velocity. 

4.2. Detonation in Fuel-Ash Mixtures 

In the absence of t his viscous ignition (§14.11). or igni- 
tion by compressio n JArnett fc Livnlll994l: |P1 ewa et aLl 
l2004t iRopke et al.ll20071 ) , carbon detonation can only oc- 
cur in a mixture of hot ash and cold fuel in the distributed 
regime. But there can be no detonation within the mix- 
ture so long as a steady state subsonic flame exists, with 
burning balancing the average rate at which fuel is heated 
either by conduction or mixing. This is always true in 
the case of laminar flames where the thickness is much 
less than the critical mass for detonation, but it remains 
true for steady state flames in the distributed regime. 
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Even though the widths of the flames in Table 1 can be- 
come quite large, the turbulent eddy that sets the time 
scale for burning in A is itself subsonic. Burning occurs 
at a rate just sufficient to balance the advancement of 
the mixing and cannot become supersonic. 

This steady state is a fiction though, useful only for 
obtaining rough estimates for the size and speed of the 
mixed burning region. So long as the mixing region de- 
fined by the turbulent integral length scale contains many 
flames, fluctuations in the burning rate will average out. 
The situation remains closely analogous to the flamelet 
regime. Many flame surfaces combine to make a "flame 
brush" with fractal dimension D = 2.36. For the allowed 
range of length scales, the burning region moves at a 
speed given by the largest turbulent eddies and the in- 
dividual flame speeds, be they turbulent or laminar, are 
not important. 

The situation changes though when the entire integral 
length scale of the turbulence contains only one or a few 
flames (Table 1) The large eddies driving the mixing are 
random. Occasionally, a long time may elapse before a 
new eddy arrives. Within the region stirred by this large 
eddy, layers exist of nearly isother mal mixtures of fue l 
and ash. Such layers were seen by (|Lisewski et al.ll2000t) 
and were responsible for the "micro-explosions" observed 
in their simulations, but because of the small dimensions 
of the mixed flames studied in that paper by Lisewski et 
al., that burning never approached sonic speeds. 

4.3. Conditions for Detonation 

A necessary condition for detonation is that sustained 
burning inside a distributed flame width occur faster 
than its sound crossing time. This is a condition on the 
sonic length scale, 

^sonic = c s T nuc (T') < A < L. (18) 

where T' is some temperature < T < T as h in the iso- 
baric mixture. If T 1 in this equation were the same as the 
temperature in r nuc in eq. (|10p , one would require super- 
sonic turbulent motion in order to initiate a detonation, 
since that would imply 

L > A * (|) L. (19) 

The fallacy in this argument is that A is some approxi- 
mate length scale in a fictitious steady state, which never 
exists at any one place and time, while r SO nic can vary 
greatly depending on the instantaneous local values in 
a given flame. Occasionally the distribution of temper- 
ature inside the mixed flame is such that, including the 
effects of induction, it burns much faster than steady 
state. Table 3, calculated using a small 7-isotope net- 
work, gives the characteristics of burning on different 
time scales. Assume, as we shall find is necessary, that a 
high degree of turbulence exists such that L/ul ~ 0.01 
s (i.e., ul ~ 10 8 cm s" 1 , L ~ 10 6 cm). Mixing can go 
on without appreciable burning so long as the tempera- 
ture, T, remains cool enough that T nuc (T)zL/uL- Here 
a small margin of error is included, and the tempera- 
ture is calculated such that half the fuel would burn in 
0.02 s. This is To. 02 in Table 3 and the corresponding 
carbon mass fraction is Xi2,o.o2- The mixture and time 
scale are calculated for isobaric fuel-ash mixtures with 



an initial carbon mass fraction of either 0.5 or 0.75. Be- 
cause carbon burns to different compositions at different 
temperatures, this does not give exactly the same values 
as burning a given composition in place without mixing, 
but for the mixing that goes on before burning this is the 
correct procedure. 

These mixed plasmas are then allowed to burn with- 
out further intervention. Because of the high tempera- 
ture sensitivity of the reaction rate, most of the energy 
from burning is released after mixing during this phase 
of "inductive burning" . Towards the end, the burning 
can becomes supersonic for a region, r son ; c , given by the 
sound speed and burning time. The minimal burning 
time, hence smallest, r sonic , called here r™ mc , is eval- 
uated by two criteria. One is that 60% of the initial 
carbon has burned. This corresponds to the onset of the 
rapid rise in e.g., Fig.[T] A second value, the actual min- 
imum value of r son i c comes from numerically evaluating 
the point where the energy generation is a maximum for 
an isobaric mixture of carbon and ash starting at the 
initial temperature, density and carbon mass fraction. 
The initial temperature matters little. The first choice 
gives a lower temperature and hence larger r son i c and is 
thus more difficult to achieve. However, it gives a larger 
carbon mass fraction which decreases the necessary mass 
for detonation. This first set of numbers was plotted for 
Xi2,initiai = 0.5 in Fig. 4. 

For fuel densities below 7x 10 6 g cm -3 and X 12 = 0.50, 
the ash temperature is just too low to burn in less than 
0.02 s, no matter what the mixture. The sonic radius 
is much larger than the integral length scale and start- 
ing to be a significant fraction of the white dwarf itself. 
It seems that spontaneous detonation will be impossible 
below this density and a carbon mass fraction of 0.50. 
Somewhat lower densities can be accommodated in the 
carbon abundance is 0.75, but no lower than 5 x 10 6 g 
cm.- 3 . 

The sonic radii in Table 3 can be fit by a function 

r£*3c « 3.5 km F (J^j " , (20) 

where 0.4sFsl. An alternate derivation that assumes a 
constant sound speed (4000 - 5000 km s^ 1 ) and evalu- 
ates the nuclear time scale assuming T max oc T as h gives a 
similar result. 
The condition r™™ c sA in eq. (fT?) gives 

/ \ 0.70 

p DDT * 1.1 x 10 7 g cm- 3 (If) (21) 

/ u L \ - 55 p-0.29 / 10 km ') - 15 
V 10 8 cm sF 1 1 V L ) 

In addition, the condition that r™™ ic be less than L ~ 10 
km gives 

Pddt * 8 x 10 6 g cm- 3 F ' 22 ^ L w - 2 (22) 

To check whether the conditions derived here are 
not only necessary but sufficient, an offline calculation 
of detonation was carried out u s ing t h e Kepler code 
dWeaver, Z immerm an, k, Wooslevl [1978; Wooslc v et al.l 
I2002D and the procedure defined in lNiemever fc Wooslevl 
(1997). Arguably, this might err on the conservative side 
since it assumes a spherical distribution and the shock 



G 



wave suffers some geometric dilution moving out, but 
the procedure has been shown to give qualitative agree- 
ment with more rigorous calculations (jDursi fc Timmesl 
2006). A region of size of variable length, L\, with char- 
acteristics given in Table 3 (T max , p m ax) was embedded 
in a comparable size region, L2 where the temperature 
declined gradually to the cold fuel value. A sample set up 
was shown in Fig. [3j It is important that the isothermal 
region is surrounded by a region where the temperature 
changes gradually so that a supersonic phase velocity 
can develop. This gradual change in temperature seems 
reasonable given the turbulent mixing that occurs on all 
length scales. 

The results for some trial detonation are given in Ta- 
ble 4 and Fig. g] It seems that homogeneously mixed 
regions of size greater than a few times r™™ c are capa- 
ble of igniting detonations. The condition r™™ ic « A 
is not only necessary for detonation, but sufficient. The 
fact that the size is somewhat larger than r son j c is not 
alarming because ignition may involve several adjacent 
mixed layers. 

4.4. Dependence on the carbon abundance 

The carbon abundance affects the transition in several 
ways, sometimes subtly. One might expect that, due to 
its larger energy release, burning a more "incendiary" 
mixture of carbon and oxygen, one with a larger carbon 
mass fraction, would make it somehow more likely to 
detonate. In fact, eq. (f22|) shows the opposite behavior, 
a weak reciprocal dependence of the detonation density 
upon the local carbon abundance. This is because the 
burning of a carbon-rich composition produces a hotter 
ash, and the fuel ash mixture maintains a thin burning 
width until a lower density. This is unavoidable. 

Unfortunate ly, this is the oppo site behavior to that 
postulated by lUmeda et al.l (|1999h in an attempt to ex- 
plain the preponde rance of bright Type la supernovae i n 
late type galaxies ((Branch et alj|199« lFirippeno~lll989f ). 
The transition densities computed here are also lower and 
the dependence on density, comparatively weak. 

However, there is another possibility. Equation (|22[) 
is a necessary, but not sufficient condition for detona- 
tion. To detonate, one also needs a critical mass. Table 4 
shows that when the homogeneously mixed region (Table 
3) is two to three times r™™ c detonation usually happens. 
Only one model is this paper strictly satisfies the criteria 
Lt ~ 3r™ n n ic < A for u L <10 8 cm s" 1 at 10 km and that is 
the model with X i2 = 0.75 and pf uc i = 6 x 10 6 g cm~ 3 . 
Turning the density down to find a solution where this 
works for X12 = 0.5 requires mixed regions that are either 
considerably larger than the integral length scale, or so 
cold that they do not detonate at all. Certainly the accu- 
racy of this study is not enough to rule out models with 
X12 = 0.5, which do satisfy this condition within a factor 
of about 2 (X12 = 0.5, p rm juei = 10 7 g cm~ 3 , for exam- 
ple, but a possibility is that detonation requires a certain 
minimum value of carbon mass fraction. This would have 
major ramifications. Models with less than critical X\2 
would have to be pure deflagrations and would have dis- 
tinctly different properties. A lot more study is needed 
before this possibility is taken too seriously. 

Alternatively, there may be other explanations for the 
preponderance of bright Type la supernova in late type 
galaxies. Perhaps, the carbon-oxygen ratio in the outer 



layers of the white dwarf at the time it explodes corre- 
lates differently with stellar evolution than Umeda et al. 
assume. In this regard it is noteworthy that more mas- 
sive white dwarfs, which are derived from larger stars, 
do have lower carbon abundan ces in their interiors when 
they explode (Fig. 6 and 12 of lUmeda et~aT1ll999f l. The 
carbon abundance, and other factors, may also influence 
the ignition conditions which are known to have a major 
effect on the supernova brightness. 

5. CONCLUSIONS 

Using simple scaling arguments, the steady state width 
and speed of a turbulent flame in the distributed regime 
have been derived as a function of turbulent energy, den- 
sity, and carbon mass fraction. This speed is much faster 
than the laminar speed at the same density, and displays 
a sensitive reciprocal dependence on the density. All such 
steady state flames are subsonic so long as the turbulence 
driving them remains subsonic. The bulk propagation of 
the burning in the distributed regime, prior to any deto- 
nation, will still be governed by the motion of the large 
turbulent eddies, not the speed of the individual flames. 

As the density nears 1 x 10 7 g cm~ 3 , however, the 
width of the flame approaches the integral length scale 
and the unsteady nature of the burning becomes impor- 
tant. For sufficiently int ense turbulence, deto nation be- 
comes possible (see also iLisewski et al.l [2000) . In order 
for detonation to occur, several criteria must be satisfied 
simultaneously. 

First, the carbon and oxygen burning flames must sep- 
arate spatially. If both fuels burn simultaneously, the 
resulting ash is too hot, and the fuel-ash mixture, too 
combustible. This makes the average width of the mixed 
burning region narrow and prevents the detonation due 
to flame broadening described here. 

Second, the speed of the largest turbulent eddies must 
approach the sound speed. That is, the maximum speed 
below which the assumption of isotropic Kolmogorov tur- 
bulence is valid must not be too subsonic. It is this 
speed that sets the characteristic scale of the problem. 
Some additional increase in flame speed may achieved 
by unsteady burning, but the amplification is not likely 
to be large enough to bridge orders of m agnitude. It i s 
the large turbulent energies reported bv iRopkel (|2007f ) . 
~ 20% sonic at the integral length scale, that makes 
things work. 

Third, the size of the mixture of fuel and ash must be- 
come as large as the largest eddies. That is, the flame 
width must approach the integral length scale. For white 
dwarfs this is some multiple of 10 km. The most likely 
detonation site will be at the merger of several such 
flames. 

These conditions may occasionally all be satisfied at 
the low densities encountered by the flame as it moves 
to the surface layers and the white dwarf begins to come 
apart. However, just reaching a particular density is not 
adequate. Nor is it sufficient simply to move into a region 
where the burning is distributed (e.g. , the Gibson length 
is sm aller than the flame thickness Ropkc fc Niemeverl 
2007). Most of the models in Table 1 are in the dis- 
tributed regime, but only those where the flame width 
A approaches the integral length scale, L, can detonate. 
Detonation thus requires the right combination of tur- 
bulent energy and density and such conditions are rare. 
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Table 1. Turbulent flame properties 



■v 


Pfuel 


ul at 10 km 


" 1 " 

1 ash 


Pash 


A 






(10 g cm °) 


(10 cm s J 


(10 9 K) 


<A[\1 — 3 

(10 g cm ° 


(cm) 


(cm s J 


0.5 


0.6 


I 


2.45 


0.284 


IJ..J.OIDIJ 


1.0(7) 


0.5 


0.(5 


3 


2.45 


0.284 


[6.15(6)] 


\5 5(7)] 


0.5 


0.8 


3 


2.64 


0.400 


5.19(5) 


2.4(7) 


0.5 


0.8 


10 


2.64 


0.400 


[3.16(6)] 


\\ 5(8)1 


0.5 




I 


2.79 


0.523 


1.55(4) 


2.5(6) 


0.5 


1 


3 


2.79 


0.523 


8.30(4) 


1.3(7) 


0.5 


1 


10 


2.79 


0.523 


4.92(5) 


7.9(7) 


0.5 


2 


1 


3.30 


1.18 


6.82(1) 


4.1(5) 


0.5 


2 


3 


3.30 


1.18 


3.54(2) 


2.1(6) 


0.5 


2 


10 


3.30 


1.18 


2.16(3) 


1.3(7) 


0.5 


2.3 


1 


3.40 


1.39 


2.43(1) 


2.9(5) 


0.5 


3 


1 


3.63 


1.89 


3.62(0) 


1.5(5) 


0.5 


3 


10 


3.63 


1.89 


1.15(2) 


4.9(6) 


0.75 


0.(5 


1 


2.68 


0.23 


3.85(4) 


3.4(6) 


0.75 


0.(5 


10 


2.68 


0.23 


1.22(6) 


1.1(8) 


0.75 


0.8 


1 


2.89 


0.33 


2.99(3) 


1.4(6) 


0.75 


0.8 


10 


2.89 


0.33 


9.44(4) 


4.6(7) 


0.75 


1 


1 


3.07 


0.43 


4.38(2) 


7.6(5) 


0.75 


1 


10 


3.07 


0.43 


1.38(4) 


2.4(7) 


0.75 


2 


1 


3.69 


1.00 


1.60(0) 


1.2(5) 


0.75 


2 


10 


3.69 


1.00 


5.07(1) 


3.7(6) 


0.75 


3 


10 


4.09 


1.62 


2.47(0) 


1.3(6) 



Whether they are sufficiently rare as sometimes not to 
happen is beyond the scope of this paper. 

ft is interesting though that the range of detonation 
densities derived here, 0.5 — 1.5 x 10 7 g cm" 3 , is what 
has been invoked for some time in order to achieve good 
agreement with nucleosynthesis, spectra, and light curves 
in artificially param etrized descriptions of the explosion 
(|Hofhch et al.lll995h . 

Many of the approximations in this paper and its con- 
clusions warrant careful checking by numerical experi- 



ments that this paper will hopefully motivate. 
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Fig. 1. — Variables in an idealized steady state flame in the distributed regime for a density 1 X 10 7 g cm -3 and turbulent speed 10 7 
cm s _1 on a length scale of 10 km. For other turbulent energies, the x-axis would be multiplied by (n^/10 7 ) 3 / 2 where uj, is the turbulent 
speed on the integral length scale (10 6 cm) in cm s~ 1 . Variables have been renormalized for plotting. The temperature has been divided 
by 3 X 10 9 K; the density, by 1.0 X 10 7 g cm~ 3 , and the carbon consumption rate, e, by 100 s — 1 . To the left of the vertical dashed line 
at the origin is unburned fuel. To the right of the vertical dashed line at 0.154 km is ash. For the assumed turbulent speeds, the whole 
distribution is moving to the left with a steady state speed of 25 km s — 1 (see Table 1). 



Table 2. Ash temperature 



p 7 X12 = 0.50 X w = 0.75 p 7 X12 = 0.50 X i2 = 0.75 



0.6 


2.45 


2.68 


2.0 


3.30 


3.69 


0.8 


2.64 


2.89 


2.2 


3.37 


3.78 


1.0 


2.79 


3.07 


2.4 


3.44 


3.86 


1.2 


2.92 


3.22 


2.6 


3.51 


3.94 


1.1 


3.03 


3.36 


2.8 


3.57 


4.02 


1.6 


3.13 


3.48 


3.0 


3.63 


4.09 


1.8 


3.22 


3.59 


3.2 


3.68 


4.16 



9 



1 7 1" 




density (gm cm 3 ) 



Fig. 2. — For a carbon mass fraction in the fuel of 50%, the figure shows the relative sizes of the sonic radius (see text) and the turbulent 
flame widths for two values of turbulent energy characterized by velocities on a scale of 10 km of 10 7 cm s — 1 (A7) and 10 8 cm s — 1 (As) - 
see also Table 1. For densities much greater than 10 7 g cm -3 , despite being in the distributed burning regime, the flame will be too thin to 
contain a region capable of running away supersonically. For lower density, however, and large turbulent energies, detonation is possible. 
No solutions smaller than 100 km could be found for r son j c for densities below 7 X 10 6 g cm -3 . For the smaller turbulent energy shown, 
detonation is unlikely. 
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Fig. 3. — Assumed conditions set up by isobaric mixing in a fuel-ash mixture with a fuel density of 1 X 10 7 g cm 3 and a carbon mass 
fraction of 0.50 (see Table 3). The nearly isothermal region in the inner 10 km has been heating from burning for about 0.01 s following 
mixing that produced a temperature, T^ax mix ~ 1-81 X 10 9 K and is now poised to runaway.The temperature in the isothermal region is 
2.20 X 10 9 K, just slightly less than the 2.26 X 10 9 K given in Table 3. The central density is 7 X 10 6 g cm 3 and the carbon mass fraction 
is 0.20. For purposes of plotting the temperature has been divided by 3 X 10 9 K, the density, by 1 X 10 7 g cm -3 and the pressure, by 
3 X 10 24 dyne cm" 2 . 
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Fig. 4. — The evolution of the conditions shown in Fig. \3\ at a time 0.68 ms later. The carbon in the isothermal region has mostly 
burned away and a detonation is forming in the surrounding temperature and carbon gradient. This detonation was followed in a separate 
calculation for over 100 km and rapidly grew in strength to consume both oxygen and carbon. The variables here are scaled by the same 
amounts as in Fig. [3] The velocity, u, is divided by 4000 km s _1 



Table 3. Conditions for supersonic burning 



Xl2 


p7 


Xi2,0.02 


To. 02 


P0.02 




T m ax 


Pmax 


T"nuc min 


„min 
sonic 




(10 7 K cm~ 3 ) 




(10 9 K) 


(10 7 g cm~ 3 




(10 9 K) 


(10 7 g cm~ 3 


(sec) 


fern) 


0.5 


0.8 


0.26 


2.09 


0.53 


0.20 


2.22 


0.50 


1.7(-3) 


7.5(5) 


0.5 


1.0 


0.36 


1.77 


0.78 


0.20 


2.26 


0.66 


7.7(-4) 


3.5(5) 


u. o 


1 ^ 


n zin 


1 fi 1 ^ 


1 97 


u.zu 


2 43 




1 lf-4) 




0.5 


2.0 


0.41 


1.60 


1.76 


0.20 


2.56 


1.46 


2.9(-5) 


1.4(4) 


0.5 


2.5 


0.43 


1.56 


2.26 


0.20 


2.67 


1.87 


1.0(-5) 


5.4(3) 


0.5 


3.0 


0.43 


1.54 


2.75 


0.20 


2.77 


2.12 


4.8(-6) 


2.5(3) 


0.5 


0.8 


0.26 


2.09 


0.53 


0.13 


2.36 


0.47 


1.2(-3) 


5.7(5) 


0.5 


1.0 


0.36 


1.77 


0.78 


0.12 


2.46 


0.61 


4.7(-4) 


2.2(5) 


0.5 


1.5 


0.40 


1.65 


1.27 


0.11 


2.70 


0.97 


6.0(-5) 


3.0(4) 


u.o 


9 n 




1 fin 




u. ±u 


Z.oo 


l.OO 


1.0[-0) 




0.5 


2.5 


0.43 


1.56 


2.26 


0.10 


3.02 


1.73 


5.1(-6) 


2.7(3) 


0.5 


3.0 


0.43 


1.54 


2.75 


0.098 


3.14 


2.04 


2.4(-6) 


1.3(3) 


0.75 


0.6 


0.43 


2.17 


0.34 


0.30 


2.36 


0.30 


6.9(-4) 


3.3(5) 


0.75 


0.8 


0.61 


1.68 


0.61 


0.30 


2.47 


0.44 


1.8(-4) 


8.6(4) 


0.75 


1.0 


0.63 


1.62 


0.81 


0.30 


2.60 


0.58 


5.2(-5) 


2.6(4) 


0.75 


1.5 


0.66 


1.55 


1.30 


0.30 


2.83 


0.93 


6.4(-6) 


3.3(3) 


0.75 


2.0 


0.67 


1.52 


1.78 


0.30 


3.01 


1.29 


1.5(-6) 


8.2(2) 


0.75 


2.5 


0.68 


1.49 


2.27 


0.30 


3.15 


1.67 


5.e(-7) 


3.0(2) 


0.75 


3.0 


0.69 


1.47 


2.77 


0.30 


3.26 


2.06 


2.4(-7) 


1.3(2) 


0.75 


0.6 


0.43 


2.17 


0.34 


0.24 


2.43 


0.29 


6.3(-4) 


3.1(5) 


0.75 


0.8 


0.61 


1.68 


0.61 


0.21 


2.60 


0.41 


1.4(-4) 


7.1(4) 


0.75 


1.0 


0.63 


1.62 


0.81 


0.21 


2.75 


0.53 


4.2(-5) 


2.1(4) 


0.75 


1.5 


0.66 


1.56 


1.30 


0.19 


3.04 


0.85 


4.7(-6) 


2.5(3) 


0.75 


2.0 


0.67 


1.51 


1.79 


0.19 


3.27 


1.19 


l.l(-6) 


5.8(2) 


0.75 


2.5 


0.68 


1.49 


2.27 


0.18 


3.45 


1.53 


3.6(-7) 


2.0(2) 


0.75 


3.0 


0.68 


1.46 


2.77 


0.17 


3.60 


1.89 


1.5(-7) 


8.5(1) 



Table 4. Detonations 



Xl2 


P7 


Xl2,l 


Ti 


Pi 


Li 


L 2 


detonates? 


0.5 


0.8 


0.20 


2.22 


0.50 


1.4(6) 


1.4(6) 


no 


0.5 


0.8 


0.20 


2.22 


0.50 


2.8(6) 


2.8(6) 


yes 


0.5 


1.0 


0.20 


2.26 


0.66 


5.2(5) 


4.4(5) 


no 


0.5 


1.0 


0.20 


2.26 


0.66 


1.1(6) 


9.4(5) 


yes 


0.75 


0.6 


0.30 


2.36 


0.30 


5.5(5) 


5.0(5) 


no 


0.75 


0.6 


0.30 


2.36 


0.30 


1.1(6) 


1.0(6) 


yes 


0.75 


0.8 


0.30 


2.47 


0.44 


1.0(5) 


7.8(4) 


no 


0.75 


0.8 


0.30 


2.47 


0.44 


2.2(5) 


1.6(5) 


yes 


0.75 


1.0 


0.30 


2.60 


0.58 


2.8(4) 


2.2(4) 


no 


0.75 


1.0 


0.30 


2.60 


0.58 


6.1(4) 


4.8(4) 


yes 



